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Abstract 

Two results are presented for reduced Yang-Milis integrals with difFerent symmetry 
groups and dimensions: the first is a compact integral representation in terms of the 
relevani variables of the integral, the second is a method to analytically evaluate the 
integrals in cases of low order. This is exhibited by evaluating a Yang-Mills integral 
over real symmetric matrices of order 3. 

1 Introduction 

In recent years the complete reduction of the D-dimensional SU{N) Yang-Mills theory 
was studied by many authors. The reduced supersymmetric partition function is an 
ordinary multiple integral of the form 

Z^y^rf^M^^e-^'^'*) (1.1) 

where the Euclidean action is 

S{X,^) = -Tr[X^,X.][X^,X.] -Tr*r^[X^, *] (1.2) 

and the matrix-valued gauge potentials X^ and their fermionic superpartners have 
values in the SU (N) Lie algebra 



[T^^r'] =ir'"T, , TTT''T'' = 5ai (1.3) 

There are several motivations to investigate these integrals both in the maximally 
supersymmetric models, as well as in cases with less supersymmetry Q We 



1 



refer to some recent papers |^ for introduction to the subject, review of recent 
results and further references. The interest in understanding these integrals is such 
that the study of simpler models without fermionic partners and for a variety of groups 
has been considered very useful, also by means of extensive numerical investigations 



|8[ |9[ |10| Unfortunately, a general rigorous method for computing eq. dl-lD 

is still unknown. The development of techniques to understand these integrals for 
various values of A'^ and D is very desirable, as well as techniques which allow the 
analytic evaluation of the integrals in cases of small order of the matrices. 
In this letter we concentrate on the simpler, pure bosonic integrals of the form 



Z= / ^dX„e^l^-^"'[^-^-l (1.4) 



The existence conditions and the convergence properties of (1.1) and ( |l.4| ) (and of the 
corresponding correlation functions) have been analytically established in the most 
general case and for any compact semi-simple gauge group recently in As it is 
well known, these models, besides the global unitary invariance for similarity transfor- 
mation of the matrices, have a rotational SO{D) vector invariance. The latter is the 
crucial ingredient for our two results : the first is a compact representation of eq. ( |l.4[ ) 
in terms of its natural variables. It holds for generic values of D and any group. The 
second result is a technique to analytically evaluate the integrals in the case of small- 
order matrices. We exploit it in the evaluation of a 5D integral corresponding to the 
Y-M integral over 3x3 real symmetric matrices. It seems likely that both results 
should be helpful also for the supersymmetric integrals. 

It is useful to parametrize the random matrices with their independent entries. For 
instance we shall consider the ensemble of real traceless symmetric matrices of order 
3 : 



a" 


u" 


v'" 






z" 




z" 





K" = i u^" -a" -6" 2" ; n = l...D (1.5) 
V z" &^ / 

The action S{X) is easily evaluated as a polynomial in the scalar products of the five 
independent D-dimensional vectors. The ensemble has the unitary invariance S'0(3) 
and we evaluate, for _D > 6 

ry I jD ,Dr,D ,D jD -S(X„) 

^so(3) = / " o^" O" ud v a z e ^ 



(z?-4)(i3-6)r(Mz_i)r(f-i; 



(1.6) 



This should not be confused with the ensemble of real anti-symmetric matrices, also 
having the same SO(3) unitary invariance. The latter ensemble is parametrized by 
three D-dimensional vectors and corresponds to the only known integral which was 
exactly evaluated in every dimension both in the bosonic and supersymmetric case Q 
|5| |l^] . Further evaluations we re g iven for supersymmetric integrals in fixed dimension 
by the deformation technique ||l^, reproducing the conjecture of ||l^ and afterwards 
numerically confirmed for small gauge groups in 0. 
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2 A compact natural representation 



Let us consider an ensemble of matrices X,,, fi = 1, . . . , D, parametrized by a number 
n of D-dimensional real vectors. For instance, in the familiar ensemble of Hermitian 
traceless matrices of order A'^, invariant under the global SU (N) group, there are 
n = N"^ — 1 real vectors. The action S — — Tr is a homogeneous polynomial 

in the components of the vectors, and it only depends on the n(n + l) /2 scalar products 
of the n vectors, as it foUows from the vector SO{D) invariance of the model. The 
bosonic partition function is a multiple integral over nD variables, with an integrand 
which is a function of n{n+l) /2 variables. If n < D the vectors are generically linearly 
independent and it seems desirable to express the partition function as an integral over 
the natural variables, that is the scalar products. This is expressed by the equation 

Z ^ [Y[dv, F{B) = f Y[db,, F{B) In,D(B) , 

bij = Vi ■ Vj = {vi\vj) ^ Ulj Zij , li = \/ (viY (2.1) 

Here F{B) stands for the usual Boltzmann weight (possibly multiplied by a Pfafhan 
or a determinant resulting from the integration over the fermionic fields), the region 
of integration R is described later in this section and In.oiB) is the joint probability 
density for a set of n vectorf 
dimensional Euclidean space. 
Proposition 1 : 

Iu.d{B) = 




where no-i = '^ ^{0/2) ' det„(B) = (/i^a . . . ^n) det„(Z) (2.2) 
Proposition 1 may be proved from^ 

Proposition 2 : for a given set of n— 1 linearly independent vectors {vk}, k = 1, .., n— 1, 
in a real D- dimensional space, n < D, 



dVn S \l — ] 5 \ Zln " • • • Zn^l,n ; T 

no-n D det4Z)(^— ^)/^ ^ no-n „+1 /W \ det„(i3)(^— 1)/^ 
2 " det„-i(Z)(o-")/2 2 " Ili'-'J det„-i(B)(-D-")/2 



(2.3) 



differcnt proof of eq. (2.2) may also be given by writing the integral representation for 
the delta functions and evaluating the integral of a random real symmetric matrix in external 
field. 
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where det„_i(i3) is the leading principal minor obtained by deleting the last row and 
last column of B. 

We outline here a proof for Proposition 2. The vector v„ may be decomposed on the 
basis of the n — 1 external vectors and a D — (n — l)-dimensional vector Vn± in the 
orthogonal subspace : 



duji 



(2.4) 



The vector v„± only occurs in the first deha function 5(1 — v^/l^), and its integral is 
trivial 



V„± (5 1 — 



dpp 



5 1- 



~ f' n. 



n-l 

E 



LU j LU k bjk 



(2.5) 



The remaining n— 1 delta functions fix the lu j as the solution of the system Ej-i ^'-i 
Idn Zin — bin, i = 1, .., n — 1, that is LUj = '^^[B~^]jsbsn. Then 



2 n — i 



UJjUJk bjjc — bnn - 



n-l 

E' 



-h h [R-il /, _ detn(^) 



Finally Proposition 2 is proved after using 

n — l n — l / 1 



fc=i 



E,— 1 



In 



('n 



/1/2. -L-i det„_i(Z) 



(2.6) 



(2.7) 



Proposition 1 foUows now from Proposition 2 since the joint density J„,_D(i3), from its 
definition, may be written as the product of the integrals ik- 

The cosines Zij are bound by — 1 < Zij < 1 but the derivation of eq. (2.2) implies 
the stricter bound that detfe(Z) > O for each k, k — 1, .., n, which is necessary and suf- 
ficient condition for the real symmetric matrix Z to have all the eigenvalues positive. 
Eqs.(|^), (U) provide a very simple reduction of the partition function to the smaller 
number of relevant variables which has general validity for different matrix ensembles 
and dimensions. It may be useful for numerical integration and possibly for more gen- 
eral purposes when averaging over random vectorsn The partition function expressed 

■^See for instanccLchapter 21 of ref . [p^ , where one studics averages over randomorthogonal 
vectors. Our eq. (2.2) for n = 2, Z12 = O , reproduces the result eq. (21.1.11) of llSl. 
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as multiple integral over the scalar products is closely related to the Ingham-Siegel 
integral |l^. In lact the well known Yang-Mills integral for the SU{2) invariant en- 
semble is a trivial case. After parametrization ol the hermitian traceless matrices of 
order 2 with 3 real entries, one obtains 



3 

„. „-Ss 
"J 



ZsU(2) = /"j^'^^'j 
J = l 

^3{D-l)/2 p 

dB e"''s[/(2) [det(B)]^ 



where 



'S's(7{2) = 16^&ii&22 + &11&33 + &22&33 - &12 - bi3 - &23^ = 16det(i3) tr(_B ^) . 

We change integration matrix variables dj = (det _B)[_B~^]ij , dB = (detC)~^dC, 
detB = (detC)^''^ 

dB e-^sc.(2, [det(B)l(°-*'/^ = 2"^° / dC (det e"*"-^ = 

3/2 (D\^(D-1\ (D-2\ , , 



23D 



Since for the 517(2) case in 73 = 10 the Pfafhan arising from the fermion integration 
is proportional to (det B)'', also the supersymmetric SU{2) integral is easily obtained, 
reproducing the result of Q p^ . 



3 Reduction by iterated projections 

A useful technique which avoids the inconvenience of the Jacobian and positivity 
requirements is the Reduction by iterated projections. The price is a more complicated 
expression of the action. To put it on general ground, let us consider the integral in 
n real vector variables of a function which only depends on the scalar products in D 
dimensions 



/= / d'^Xld°X2 ■ ■ ■d'^TnFi {xi\xj)D;i < j ] (3.1) 



Proposition 3 
Case D > n: 



I = ■ ■ - f^D-n / dxiXj^^... / dx„x°~''^ / JJ^d^ij F(. . .) (3.2) 



{xi\xj)n = xiiXij + . . . + Xi-i^iXi-ij + XiXij, i < j (3.3) 
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Case D < n: 



1 = 0., 



■ Oo 



dxr 



Yl dxij F{. 



(3.4) 



i<j 



{Xi\Xj)D = X\iXlj + . . . + XDiXDj, D<i<j ■ 

{xi\xj)D = xuxij + . . . + XiXij, i < j, i < D (3.5) 

Proof; The rotational invariance in the space of all vectors allows us to choose 
the vector xi as the polar direction. Any other vector Xi is then described by its 
component parallel to xi, denoted as xu Uving in R, and a vector orthogonal to it, 
which (with abuse of language) we again denote Xi, Uving in the space R^~^. Thus 
we have the first sweep 

P CO poc f 

I = Od-1 / dxix^~^ j dxi2...dxin / dx2~^ ■ ■ ■ dx^~^ F{. . .) (3.6) 

JO J-OD J 

where the scalar products in F are easily evaluated: 

{xi\xi)d=xI, {xi\xi)D = xixii, {xi\xj)D = xuxij + {xi\xj)D~i, l<i<i 

Next we choose the vector X2 as the polar direction in R^~^ , and introduce accordingly 
new real variables X2i and vectors Xi in R^~^, i = 3, . . .n. The integral is now: 

/•oo p oo poo 

I = ftD-lOD-2 / dxixf'^ / dX2X2~^ / dxi2 ■ ■ ■ dxindX23 ■ ■ ■ dx2n 
JO Jo J -oo 

d°-\3...d°-\^F{...) (3.7) 
The scalar products in R^~^ that still appear in F are evaluated: 

{x2\x2)d-\ = X%, {x2\xi)D-l = X2X2t, {Xi\Xj) D-\ = X2iX2j + {Xi\Xj) D-2, 2<i<j 

The process is iterated and we have two situations: D > n, D < n. In the first case we 
can iterate the process n times, to obtain an integral in n positive variables xi ... Xn 
and in(n — 1) real variables Xij, i < j. When D < n the process can be iterated only 
D times; we obtain D positive variables Xi ... xd, and ^D{2n — D — 1) real variables 
Xtj, i < j, i = 1 . . . D. 

4 Y-M integral with 3x3 real symmetric matri- 
ces 

The 3x3 traceless r eal s ymmetric matrices of the Yang-Mills integral are parametrized 



as indicated in eq. (1.5). The partition function 

Zso(3)^ I d°ad°hd°ud°vd°ze'^^'> (4.1) 



6 



contains the foUowing expression for the action, in terms of scalar products like 



So = —4:[a^u^ — (a|u)^] — [a^ — {a\z)''] — [a^v^ — {a\v)'^] — [b^u^ — (b\u)^ 
-4[6^2^ - {b\zf] - [6^t;^ - {b\vf] - [uz^ - {u\zf] - [v^z^ - {v\zf] 

-[u\'' - {u\vf] - 4[M''(a|6) - {u\a){u\b)] - ^z\a\b) - iz\a){z\b)] 
+2[v\a\b) - {v\a){v\b)] + 6{a\z){u\v) - 6{a\v){u\z) 
+6{b\u){v\z) - 6{b\v){u\z) 



(4.2) 



The huge invariance of the action can be exploited by reducing to a smaller number 
of eflfective variables. We have been able to compute Zso(3) by using the technique of 
successive projections, described in the previous Section. 
For D > 5 the reduction yields a 15-dimensional integral: 

poo poo poo foc poo 

Zso(3) = 0,D-i ■ ■ - i^D-s / dvv'^~^ / dzz'^~^ j duu^~'^ / daa^~'^ / dbb'^~^ 
Jo Jo Jo Jo Jo 

dzv duv dttv dbv duz daz dbz duu dbu dba e'^^"'"'^^"'"'^^"'"'^'' (4.3) 



J —c 



where the action in the reduced variables has been decomposed in terms with the 
foUowing meaning. Si includes all terms in the variables o, ba and b, that terminated 
the reduction procedure and can be integrated immediately: 



Si = ~4a^[v^ + iz'' +zl) + 4:{u +ul + ul)] 

-4(fe2 + bl)[v^ + 4(2^ + zl) + («' + ul+ ul)] 
+8aba[v'^ — 2{z^ + z^ +ul + ul)] 



(4.4) 



The integration yields a factor that does not depend on the 6 variables ax, bx, 
where x = u,v, z. These variables appear in the action quadratically and linearly. 
The quadratic terms are collected in S2 and the linear terms in S3, in a form suitable 
for a Gaussian integration: 



S2 = -4X*MX 



/ Uv + 2bv \ 
2ov + bv 



where 



X = 



M : 



ZZy 

O 
O 

\ o 



o 

+ u1 
O 

O 

—UUv 



dz + 2bz 
2az + bz 

ttn + 2bu 

\ 2au +bu/ 
O 



^zi 

-v'" 

o 
o 



5*3 = 8vX*y 



Y 



-UzUv 



(4.5) 



2 I 2 , 2 
v + u + w„ 

o 

-UUz 





— ZUz 


\ 




-ZUz 








S'iijj; Zl} ZUu 








UzZy ~\~ ^ZUy 








2uzv 






v 


— UZv 


J 





















— UUv 



















— UUz 




+ + 4 







+ul + ul/ 
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The remaining part 54 of the action contains terms in the reduced variables originating 
from the off diagonal matrix elements, namely v, z, «„, u, Uy and Uz'- 

S4 = —AZ'^{V'^ + M^) — 4:Ul{V^ + Zy) — 4:U^{V'^ + Z'^ + Zl) + iZUzUyZy 

(4.6) 

As we described, the Zso(z) integral can be done in independent steps. We have a 
factor coming from the triple integral: 



Z' oc /* 00 Z' 00 

Jo Jo J — 00 



— (12)"°+^ ^( 2^) r (^2'') 



(4.7) 

At the same time one can do the Gaussian integral 



I2 ~ I daydbvd/izdbzdaudbu e 



S2+S3 



J d = ^(detM) e 



(4.8) 



The computations for inverting the matrix M and evaluating the quadratic form S5 = 
4:v'^Y*M~^Y were done with the aid of Mathematica: 

S, = 4{v' + zl){u^ + ul) + ^z^v" + ul) - Szzyuyuz - 4^^ - 4 ^^^'fl 2 - 

z^ + zi + ut+ ui 

J , T.r 4 2 2/ 2, 2\/ 2 , 2 I 2\ 

det M = v z u {z + Zy)(u + w„ + Wj) 
In summing the expressions S4, and Ss, to obtain the residual action, many terms 
cancel: 

c. , Q ,1^2 2 ^„,2^2„ 2 + ul ul + Z^ + zl , , 

64 + 65 = —4« U — 4V Z U „ J— 2\ l^-^) 



We arrive at the foUowing stago of computation of the fuU integral: 

Zso(3) = j77^/'l2-^r - 2) F - ^) «o-i . . . n^-s /3 (4.10) 
where we still have to evaluate 

poo poo poo POO 

/3 = / dvv'^~^ / dzz'^~^ / duu^~'^ / dzv duv duz e^'^'^^^ 

Jo Jo Jo J -00 

[v'^{u' +ul +ul+ + zl) + {z^ + zl){u^ +ul + W^)]~ 2 

r/ 2 . 2 . 2n/ 2 . 2\i-l/2 

[(u +Uy){Z + Zy)\ 

Since the variables w„ and Mz always appear in the combination ul + ul = p^, we 
introduce polar variables obtaining a factor 2-k from angle integration. It is also 
convenient to change from the variable v to the variable x: 



u-^ + p-' + z^ + zi 

(w2+p2)(^2+^2) 
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-'o -'o 

pdp dz^{u^ + + + zl)-^+'' (4.11) 



The double integral in p and z^ is elementary, and poses the restriction D > 5 for 
convergence, 



dz^ / pdp{u +p'' + z'' + zl)-^+^ = {z^ + u)-^ + ^B ( 

-oo J O ^ 



D-S 3 
2 ' 2 



The remaining integrals are easily done with the substitution ti = (1 + x^)~^ ^ r cos O 
and 2: = (1 + x'^)~^^'^r smO . We obtain, with the further restriction D > 6 for conver- 
gence, 

2^B(:^,^)^°°dxx°-^(l + :r^)-|^+i jT'^'decos^-^esin^-^^ 



/3 

7 U— L — 4r' sin" W cos" W < tj/^r»— ■ir 



D-l -4r^sin^ O cos^ ^ — A 5/2o-4-D 2 2 ^) 



i " ' ~ r(|D-|) 

(4.12) 



The final result , for D > 6 is: 

^ P-4)(i?^-6)r(Mz:|)r(f-i) ^'-''^ 

In the case D = 4, one should truncate the reduction by projections, as described in 
the previous section (since the 5 vectors are not linearly independent). This produces 
a set of 14 integration variables and an action which, in simple terms, are obtained 
from eq. (4.3) by suppressing the integral in b and all terms containing 6 in 5*1. We 
checked that also in this case the integrations lead to a divergent result. 
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